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The mid-infrared spectrum of the transiting exoplanet 

HD 209458b 

M. R. Swain^, J. Bouwman^, R. L. Akeson^, S. Lawler^, & C. A. Beichman^ 

ABSTRACT 

We report the spectroscopic detection of mid-infrared emission from the tran- 
siting exoplanet HD 209458b. Using archive data taken with the Spitzer/IRS 
instrument, we have determined the spectrum of HD 209458b between 7.46 and 
15.25 /xm. We have used two independent methods to determine the planet 
spectrum, one differential in wavelength and one absolute, and find the results 
are in good agreement. Over much of this spectral range, the planet spectrum 
is consistent with featureless thermal emission. Between 7.5 and 8.5 fim, we 
find evidence for an unidentified spectral feature. If this spectral modulation 
is due to absorption, it implies that the dayside vertical temperature profile of 
the planetary atmosphere is not entirely isothermal. Using the IRS data, we 
have determined the broad-band eclipse depth to be 0.00315 ±0.000315, imply- 
ing significant redistribution of heat from the dayside to the nightside. This work 
required development of improved methods for Spitzer/IRS data cahbration that 
increase the achievable absolute calibration precision and dynamic range for ob- 
servations of bright point sources. 

Subject headings: planetary systems — stars: individual (HD 209458) — meth- 
ods: data analysis — techniques: spectroscopic 



Introduction 



The Spitzer Space Telescope has revolutionized the observational characterization of 
exoplanets by detecting i nfrared emission from these obie cts; measurements have be en re- 
ported for HD 209458b jDeming et ah Il2005ah . TrES-1 Jcharbonneau et al. I ^m^ . HD 
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189733b flDeming et al. Il2006h . and v Andromeda b flHarrington et al. 1120061) . HD 209458b 



the first reported transiting exoplanet (jCharbonneau. et al. 



2000^ . is located at a distance 



Th e most recent sy s tem p aram 
(|2Qo3, 



with 



of 47 pc and has a GO stellar primary (V = 7.6 mag) 
eters for this hot Jovian exoplanet have been established by iKnutson et al. 
P = 3.52474859±0.00000038 days, Mpianet = (0.64±0.06)Mj, an eccentricity consistent with 
zero, and Rp = (1.320 ± 0.025)i?j, 10-20% larger than predicted by irradiated planet mod- 
els. The detection of infrared emission from hot Jovian exoplanets has stimulated extensive 
theoretical work on the atmospheric structure and emission of these planets. Constraining 
the model predictions for infrared emission from hot Jovian atmospheres is an important 
motivation for current observing programs. 

Spectral characterization of hot Jovian exoplanets is a high priority and is essential for 
understanding atmospheric composition and propert ies. Spectroscopic detect i on of exoplanet 



emission has proved challenging from the ground (IRichardson et al. I l2003l : iDeming et al. 



2005bl ): space-based infrared spectroscopy is particularly appealing due to the absence of 
an atmosphere, improved signal-to- noise (SNR), and instrument stability. Recently, the an- 
nouncement of a Spitze r/IRS detection of a featureless emission spectrum from HD 189733b 
JCrillmair et al. I l2007h _ and an e mission spectrum containing emission features from HD 
209458b (IRichardson et al. 1120071 ) has generated great interest. However, observations with 
the Spitzer IRS instrument are complicated by systematic errors that are large compared to 
the observable signature. Some of these systematic errors introduce wavelength-dependent 
effects; thus, careful calibration and validation is essential. In this paper we present results 
based on a new approach for calibrating the major instrument systematic effects affecting 
these observations. Using data taken from the Spitzer archive, we have determined the 
spectrum of of HD 209458b using two semi-independent methods. 



Observations 



The observations we analyzed (or iginally proposed by Richardson et al. 2007) were 
taken with the Spitzer Space Telescope (IWerner et al. II2004I ) using the Infrared Spectrograph 
(IRS; Houck et al. 2004). The data were taken on 6 July 2005 and 13 July 2005 as two 
separate Astronomical Observing Requests (AORs 14817792 and 14818048) and provide 
approximately continuous coverage of the secondary eclipse event (see Fig. [T]). The timing 
of the observations is well suited for application of the secondary eclipse technique (also 
termed "occultation spectroscopy"), in which data from portions of the orbit where light 
originates from the "star -|-planet" and "star" are subtracted to obtain the planet's emission 
(IRichardson et al. II2003I ). For both sets of observations, the IRS instrument was operated in 
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first order (7.5 to 15.2 fim) at low spectral resolution (R=60-120; SLl) with a nod executed 
at the midpoint of the observations. This observational sequence provides two completely 
independent data sets that span an interval covering the sequence: 

1. before eclipse (flux originates from star+planet), 

2. ingress (planet flux contribution changing with time), 

3. secondary eclipse (flux originates from star only), 

4. egress (planet fiux contribution changing with time), and 

5. after eclipse (fiux originates from star+planet). 

Each nod contains 140 samples with an integration time of 60 seconds each. 



To determine the orbital phase of HD 209458b we used the results by iKnutson et al. 



(120071 ) for both the period and the ephemeris. The time for each data point was determined 
using the DATE_OBS keyword in the header, which was then converted to Julian date using 
the IDL routine JDCNV.pro from the IDL astronomy library. We then converted to helio- 
centric Julian date (HJD) using the ID L routine HELIO JD.pr o (also from the astronomical 



library) for direct comparison with the iKnutson et al. I (120071 ) results. The phase was then 



estimated by phase = {{timednJrlJD — ephemeris) x mod{period)) /period. In what fol- 
lows, we will refer to the segment of the orbital phase when both star and planet are visible 
as "SP". Similarly, we refer to the segment of orbital phase when only the star is visible 
(when the planet is passing behind the star) as "S" . To determine the planet spectrum, we 
have applied the analysis described below to the spectral range of the IRS SLl module. 



3. Analysis 

The extracted fiux density time series suffers from four kinds of temporal changes (see 
Fig. [1]) that completely dominate (by a factor of ~10) th e expected signature of secondary 



eclipse fiux decrement of ~0.0025 (jPeming et al. Il2005al ). These effects are (i) a fiux offset 



between nods, (ii) a periodic flux modulation, (iii) initial flux stabilization, and (iv) mono- 
tonic flux drift within a nod. These temporal changes are not random; a scatter diagram 
shows that the flux density values are highly correlated (correlation coefficients of ~0.99). 
We find that these four major temporal flux density changes listed above are caused by (in 
order of importance) errors in telescope pointing, background subtraction, and latent charge 
accumulation. 
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Effective calibration of these systematic effects can be challenging to demonstrate. To 
test our control of the systematic errors, we developed two methods for estimating the 
exoplanet spectrum. The first method is differential and has the property that errors which 
are "common-mode" in wavelength are rejected. The second method is an absolute method 
and results in an exoplanet spectrum in units of Jy. A schematic picture of our data reduction 
method is shown in Fig. [21 central to our approach is comparing the results of the two semi- 
independent estimates of the exoplanet spectrum. Because the two methods interact with 
systematic errors differently, the comparison is useful for accessing the level of uncalibrated, 
residual systematics. In this section, we describe the initial data extraction, the major 
systematic errors, and each of our spectral extraction methods. In Section 4, we discuss the 
comparison between the differential and absolute methods for obtaining the planet spectrum. 
We then present our results and discuss the implications. We also discuss the differences 
between the methods and results of our approach and previous work. 



3.1. Data Extraction 



Our initial data extraction method is an extension of the method described by lBouwman et al. 

( 120061 ). The series of extracted images are used to define a median background image for 



each of the two nod positions. The median background image (for each nod position) is then 
subtracted from all the individual observations with the source in the other nod position; this 
generates the background subtracted images. We then identify bad pixels using a median fil- 
ter and visual inspection. The bad pixe ls are then corrected using an approach similar to the 
Nagamo-Matsuyama filtering method (INagano fc Matsuvama Ill979l). The source spectrum 
is then extracted using the method developed by Higdon et al. I ~ 2004 ) and implemented in 
the SMART data reduction package. 

The spectra were extracted using a fixed-width aperture of 6 pixels centered on the 
position of the source. The exact source position relative to the slit was determined by fitting 
a sine profile to the spectra in the dispersion direction using the collapsed and normalized 
source profile. The accuracy at which the source position can be determined is about 0.02 
pixels. This, together with the aperture width of 6 pixels, ensures that any flux variablity due 
to slight changes in the positioning of the aperture are far less than the expected planetary 
flux. 
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3.2. Systematic Errors 



Here we discuss tlie origin and chromaticity of the three significant systematic errors 
present in these data. There may be other systematic errors as well, but they, and the resid- 
uals of the errors we explicitly deal with, are smaller than the uncertainty level achieved in 
our calibration. We acknowledge that there are diff erent points of view rega.rding the calibra- 



tion o f IRS data for determining exoplanet spectra (IRichardson et al. 1120071 : iGrillmair et al. 



20071 ) and that these approaches may perform similar (but not identical) corrections to the 
data while ascribing the underlying systematics to different causes. However, ours is the 
only method that allows determining the absolute planet spectrum. 



3.2.1. Pointing Errors 
The periodic and linear drift components of the Spitzer pointing error have been doc 



umented with long-duration IRAC observations (IMorales-Calderon et al. II2006I ). Pointing 
errors cause modulation of the measured flux because telescope motion perpendicular to the 
slit axis changes the position of the stellar image with respect to the spectrometer entrance 
slit; this causes changes in the vignetting of the stellar image. Even small pointing errors 
change how the wings of the point spread function (PSF) are vignetted. Since the PSF size 
is proportional to wavelength, the measured flux changes due to pointing are wavelength 
dependent. In the absence of other effects, the measured flux density, ^(A), is 

5(A)=C(Z/,A)F(A) (1) 

where F is the "true" flux density, ( is the pointing-induced fractional flux density (C = 1 for 
no pointing error), and y is the angular error with respect to the spectrometer entrance slit 
center position in units of pixels. In principle, if ( can be determined, the effects of pointing 
error can corrected. 

We determined ({y,X), by using the spectral map observations of IRS calibrators HD 
173511 (AOR 13481216) and HR 7341 (AOR 16295168). The spectral map data consist of a 
series of pointed observations in which the star spectrum is measured on a two-dimensional 
grid (7x7 and 5x23 positions, respectively, for these AORs). For each scan perpendicular 
to the slit axis, we normalized the measured spectrum by the spectrum measured at the 
nominal slit center position, ({y) = Sy/So. Assuming the slit has a constant width, we 
combined the normalized measurements from all the slit scans. This resulted in a series of 
values at each nominal pointing position perpendicular to the slit axis {y G [yi, 2/2, 1/3, ■ • •])• 
The difference in these values at each nominal slit scan position reflects a pointing error that 
can be corrected for in an iterative process. We defined a "template" by taking the average 
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value of the points at each pointing offset position. The individual slit scan data were then 
shifted in the horizontal axis and renormalized to minimize the value of the shifted curve 
with respect to the template. After all the slit scans had been shifted and renormalized, a 
linear interpolation was done to find revised values for the template function at the nomi- 
nal pointing offset positions transverse to the slit axis. The individual slit scan data were 
then shifted and renormalized again for a best fit to the revised template function. This 
process was iterated until convergence was reached; it resulted in pointing-error-corrected, 
slit-scanned data. We determined C by fitting a cubic-spline at each wavelength through the 
shifted and renormalized slit scan measurements (see Fig. |3]top). 

While a periodic pointing error component is frequently s een in Spitzer observations, it 



is not necessarily repeatable in terms of shape or amplitude (ICarey 1120071 ). The IRS data 
we analyzed for HD 209458 show a periodic modulation of the measured flux (see Fig. [1]) 
that could be due to the Spitzer pointing error. To test the hypothesis that changes in the 
measured flux are due to pointing errors, we modelled the pointing error periodic motion in 
both the spatial and spectral axis. This leads to an elliptical motion that creates a symmetric 
profile about individual maxima and minima. The asymmetric profiles in these data require 
the addition of a harmonic term for angular velocity; when this is incorporated, the pointing 
error is given by 

e = eQ + Aesin{ujt-(t)e), (2) 

X = Xo + rrixt + AxCos[uj9(t) — (j)^], (3) 

y = yo + rriyt + AyCos[u6{t) - (py], (4) 

where t is time, x is the position parallel to the slit axis (the spatial dimension on the array), 
y is the position perpendicular to the slit axis (the spectral dimension of the array), Xo and 
yo are initial offsets, and ruy are the linear drift terms, 6 is the angular velocity, A is the 
amplitude, u is the frequency, and (p is the phase. The normalization of t and 6 is determined 
by the conditions 

t G [0, 27r], / 9dt = 2n, Ag < Oq. (5) 
Jo 

We determined the parameters for the x and 6 components of our pointing model by fitting 
to the source motion along the slit axis using the following steps for the data in each AOR: 



1. Determine position: For each measurement in the time series, we constructed the 
spatial profile at each spectral channel. These profiles are normalized by wavelength in 
the spatial axis and shifted so that they can be "stacked" coherently. A median spatial 
profile is then determined. This median spatial profile is then fit with the function 
smc^(x). The fitted position of the maximum of sinc^{x) as a function of time is used 
as the measure of telescope pointing changes in the spatial axis. 
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2. Linear fit: We fit and removed the linear component, mx, of tlie source position in 
tfie spatial axis of the data in each nod. 

3. Determine the frequency: To determine the frequency, u, of the periodic oscillation 
in each nod, we took the Fourier transform of the linearly detrended position function 
in the slit spatial axis. The normalized frequency values were the same within the 
errors, and the mean frequency was used in the remaining analysis. 

4. Characteristic profile: We folded the data, computed a median profile and local 
standard deviation, applied a 10 cr clip to remove discrepant points, and determined 
the mean profile. 

5. Determine model parameters: We determine the values for Xo, Ag, A^, (pe, 4>x by 
fitting the predicted position along the slit axis to the measured source position along 
the slit axis (see Fig. [3]); the values for these parameters are given in Table 1. 

At this point, we only need to determine the values of {yo, niy, Ay, (py) to completely describe 
the telescope pointing. Judicious selection of values for (?/o, ^y, Ay, (py) produces an estimate 
of the cross-slit position changes that (with appropriate normalization) agree remarkably well 
with the intensity time series (see Fig. [3] bottom) and successfully reproduce the asymmetric 
component in the shape of the periodic modulation. 

The excellent agreement between our simple pointing model and the observed changes 
in the measured flux confirm that, in the case of a point source, the IRS measurement of the 
flux density is affected by the position of the (stellar) image in the spectrometer entrance slit. 
The results (see Fig. |3] bottom) imply that pointing changes as small as ~ 10 milliarcseconds 
have an effect on the measured IRS flux for bright, point-like objects. Because the PSF size 
is a function of wavelength, the pointing error effect is chromatic. The asymmetry in the 
PSF wings also causes pointing errors to be asymmetric with respect to the nominal center 
of the slit (this can be seen in Fig. H]). 

Equipped with equations 1 and 4, we can now decompose the changes in the measured 
flux density into three specific kinds of pointing errors, all of which can be seen in Fig. [TJ 
Each of these pointing errors contributes a specific component of y. Note that the values of 
Ay & 0y are the same for all AOR/nod combinations, while yo & my are different for each 
AOR/nod combination. 

• initial peakup/nod error - The pointing error associated with the initial peakup or 
nod operation. The high-accuracy peakup, used for these observations, has a 1-a error 
circle radius of 0.4 arcseconds. This translates into a flux uncertainty of ~ 5 — 10%. 
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When a nod is executed, there is significant motion perpendicular to the sht axis. This 
is the reason why the median flux density differs in each AOR/nodl. The initial error 
is static and represents a constant offset described hj y = f{yo). 

• pointing drift - During an observation, there is a slow linear drift in pointing during 
each nod. The drift rate is larger at the nod2 position. The slow pointing drift rate 
ranges from 3 mas/hr to 19 mas/hr (based on a 1.85 arcsecond per pixel plate scale 
and a nod duration of 2.9026 hr - see Tab. 1). This linear pointing error is described 
by y = f{my). 

• periodic error - The Spitzer telescope has a known periodic pointing error ~ ±30 
milliarcseconds. This is the error that causes the clear periodic modulation of the flux. 
The periodic position changes are described hj y = f{Ay, (py). 



3.2.2. Background Correction 

In the mid-infrared, accurate measurement of the infrared source flux requires subtrac- 
tion of the background due to local zodiacal emission. To remove the background contribu- 
tion to the spectrum, we construct and subtract a median background image. However, this 
median image must be constructed with care as there is a systematic error in the backgound 
estimate due to leakage from the bright source. This leakage is manifested as a flux density 
offset between the background at the nodi and nod2 positions. In principle, this offset could 
be caused by structure in the background. However, inspection of IRS calibrator star data 
shows that the difference in the background between the nods is systematic in that it occurs 
for all the multiply observed IRS calibrators we checked; the effect is highly repeatable and 
is proportional to the measured source flux. IRS calibrators observed with a series of slit 
offsets show the measured source flux decreases with the slit offset from the target, and the 
background offset is proportional to the measured flux. This suggests that some of the light 
from the source contaminat es the background through t he wings of the PSF. Because the 



Spitzer PSF is asymmetric (IBayard fc Burgarolas II2004I ). the leakage differs in nod 1 and 
nod 2. 

To determine the amount of a point source contamination in the background estimate, 
we have used observations covering an interval of approximately three years for five IRS 
calibrator stars (HR 6606, HR 7341, HD 166780, HD 173511, HR 6348), together with the 
assumption of a locally uniform background. The IRS calibrators we selected were observed 
in the nominal nodi and nod2 positions for both SLl and SL2 modes. These stars were 
observed on a regular basis throughout the Spitzer operational period up to the time of 
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these observations. Each star was typically observed at least 20 times over a three year 
interval. Thus, slit precession over a period of one year is a strong test of our assumption of 
uniform background. 

We determined the source contamination in the background by subtracting two SLl 
background positions when the star is in the SLl and SL2 positions. The background source 
contribution function, BSCF, for the nodi position has the form 

BSCF. d - ^^'"'''^^"^^^ - Bnodi{SLl,nod2) - Bnodi{SL2,nodl) / RSRFnod2 \ 
Ssource{nod2) S' {nod2) - Bnod2{,S L2, Hodl) \RSRFnodi)' 

where Bnodi{.SLl, nod2) is the background at the nodi position measured when the source is 
located at the SLl, nod2 position; Bnodi{SL2,nodl) is the background at the nodi position 
measured when the source is located at the SL2, nod2 position; S'{nod2) is the measured 
source flux at the SLl, nod2 position with no background correction; and RSRF is the 
relative spectral response function at either the nodi or nod2 source position. For each 
term, the subscript denotes the position on the array where the value was measured while 
the source position at the time of the measurement is indicated by the parenthesis. Thus, 
BnodiiSL2,nodl) is the background measured at the SLl, nodi position when the source 
is located at the SL2, nodi position. Since we are calibrating SLl data, the background is 
always measured in the SLl slit. However, determining the BSCF requires using data when 
a star was observed with both the SLl and SL2 slits. The BSCFnodi is similarly defined 
except that all nodi instances become nod2 and visa versa. The corrected background flux 
density at the two SLl nod positions is then 

Sinodl) = [S'inodl) - Bnodi{SLl,nod2)] + BSCF^odi x S{nod2), (7) 

and 

S{nod2) = [S'{nod2) - Bnod2{SLl, nodi)] + BSCFnod2 x S{nodl). (8) 

The system of linear equations is then solved for the background corrected, measured source 
flux density, S*, at each nod position. This results in an estimate of the BSCF each time 
the calibrator stars were observed. We then averaged the results for all the calibrator ob- 
servations to determine a mean BSCF (see Fig. EJ. The uncertainty in the BSCF at each 
wavelength was determined by the standard deviation in the mean. 

Applying the BSCF substantially reduces the background flux density offset between 
the nodi and nod2 positions (see Fig. [5]). For wavelengths shorter than 13 ;um, the correction 
we derive is ~0.4% for nod 1 and ~0.9% for nod 2. This means that there can be ~ 0.5% 
of the source flux present in the wings of the PSF ~ 20 arcseconds away from the observed 
source position. Thus, the contamination of the background by the source is of the same 
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magnitude as the signal from the exoplanet. This correction for source contamination of 
the background may not be necessary for many observations. However, for high dynamic 
range measurements on bright point sources, neglecting the leakage of the source into the 
background estimate introduces systematic errors in the data for each nod. 



3.2.3. Latent Charge Accumulation 



In the context of the IRS instrument , latent charge accumula t ion has been reporte d 



by several authors (iGrillmair et al. 1 120071 : iRichardson et al. 1 120071 : iDeming et al. 1 120061 ) , 

and is sometimes termed "charge trapping" . Currently, the details of the semi-conductor 
physics that produce the effect are not well understood. Empirically, the responsivity of a 
pixel initially depends on the illumination. When the flux density time series is median- 
normalised (e.g. Fi(X)/ < F{X) >), the light curves at each wavelength can be "stacked" 



(IRichardson et al. 1 120071 ) . The effect of latent charge accumulation can be seen at the 



beginning of each AOR in Fig. [U the effect is characterized by a rapid initial increase in 
the measured flux density, which then approaches an equilibrium. If one excludes the first 
~ 20 points in each AOR and finds the slope of a best fit line to the data, the slope in nod2 
is greater than the slope in nodi . As Fig. [3] shows, a simple pointing model explains the 
changes in the measured flux after the first ~ 20 minutes. Note that it is possible to confuse 
the linear component of the pointing drift with the effect of latent charge accumulation 
after the first ~ 20 minutes. By explicitly modelling the pointing, our analysis breaks this 
degeneracy and allows us to separate the effects of these two systematic errors. We conclude 
that latent charge effects are negligible after the first 20 minutes. We omitted the data 
affected by latent charge accumulation from further analysis so that the effects of latent 
charge do not impact our estimate of the planet spectrum. 



3.3. Spectral Response Function 

After the initial extraction, the data were background corrected using the background 
correction discussed above. The next stage of the calibration was to derive and apply a 
spectral response function. Using the IRS calibrator i] Dor, we derived our own spectral 
response function, the SRF, for the nominal nod positions and extraction aperture. We 
selected this source for defining the spectral response function because it has the same 
brightness as HD 209458 and thus should minimize any remaining instrument residuals. 
Both HD 209458 and the r] Dor data were extracted using identical methods, and both 
incorporate identical methods for background correction. Thus, the treatment of both the 
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calibrator and source data sets is fully self-consistent. In the case of AOR 14818048, an 
additional calibration step for the SRF was required because the observations were not 
carried out at the nominal nod positions. To determine the changes in the SRF for other 
(but still relatively nearby slit positions), we used observations of the IRS calibrator stars 
HD 42525 and HR 7341, which were observed at intervals along the IRS slit. We interpolated 
between observing positions to determine how the calibrator star's spectrum changed as a 
function of slit position and used this information to renormalize the SRF derived, using rj 
Dor for the nod positions used in AOR 14818048. 

At this point, the data were ready for extraction of the spectrum. Of the three ma- 
jor systematic errors, the background had been removed (at this point in the calibration 
sequence) by explicit calibration. The effects of latent charge were removed by excluding 
the effected data from the spectral estimation. However, the pointing error remained uncor- 
rected. In what follows, two methods were used to correct the pointing error and extract 
the planet spectrum. 



3.4. Spectral Estimation 

3.4.1. Differential Method 



This approach assumes that changes in the measured flux have a wavelength-independent 
component, characterized by G{t), which can vary on a timescale of minutes, and a wavelength- 
dependent component, G{X), which is stable for a given nod but can change between 
nods. The G{X) term is removed by construction of a spectral fiat. We derived a spec- 
tral fiat for each nod by comparing the average fiux in each spectral channel to the fiux, 
-F(A), of a stellar photosphere model for HD 209458 (IKurucz I Il992l ) normalized to the 
12 fim fiux. Thus at each wavelength. A, the spectral fiat was defined as the inverse of 
[Ss{\) / FKuruczW] X [FxMrucz(12)/S'5(12)] whcrc Ss is the measured fiux observed in inter- 
val S. After normalization of each nod by the associated spectral fiat field, the data were 
as sumed to vary only i n tin ie; a more extens i ye dis cussion of this technique can be found 



m 



Bryden et al. 



(I2OO6I ) and lBeichman et al. I (120061 ). To reject the wavelength- independent 
G{t) term, we constructed a differential observable using the following method. During 
period SP (star+planet), the measured fiux, Ssp{X), can be written as 



Sspi\) = Git) [FsiX) + Fp{\)] 

This can be expanded as 



(9) 



SspiX) = Git)Fsi\) 



1 + 



FpjX) FpjX') FsjX') 
Fp(A') Fs{\') FsiX) 



(10) 



where F{X) is the true source flux, the subscripts refer to the star or planet, and A' is 
a reference wavelength selected for the comparison. We set the transit depth at A' to a 
plausible value such that {3 = {Fp{X')/ Fs{X')) << 1 and the transit depth at A, relative to 
the transit depth at A', is a = {Fp{X)/ Fp{X')). Ssp{X) can then be expressed in terms of a 
and P as 



Ssp{X) = G{t)Fs{X) 



1 + ap 



FsjX') 
Fs{X) 



During period S (star only), the measured signal, 5*5, is ^^(A) = G(t)Fs{X). The ratio 
of the two wavelengths during the SP and S periods is Rsp = Ssp{X)/ Ssp{X') and Rs = 
Ss{X)/Ss{X'). The advantage of taking the ratio is that the wavelength-independent gain 
term, G(t), drops out. Appropriate substitution, and solving for a, yields 



a=[Rsp{l + f3)-Rs]/P. (12) 

The observables are Rsp and Rs, and a is the measure of the brightness of the planet at A 
compared to the planet brightness at A'. The results in Fig. [10] reflect a value for (3 = 0.003. 
However, the results for the spectral slope are not strongly dependent on the assumed value 
for P, and we explicitly measured the echpse depth in any case (using the absolute method). 

We summarize the steps in the differential method as follows: 



1. Treat each AOR and nod combination as an independent secondary eclipse measure- 
ment with an independent calibration; this leads to four independent estimates of the 
planet spectrum. 

2. Normalize the data in a given nod by dividing the flux density by the 12 /xm value at 
each sample in the time series. 

3. Average the SP and S intervals in the time series to get the "star only" and "star+planet" 
spectra. 

4. Normalize a Kurucz model for HD 209458 flux density by the model's 12 fim prediction. 

5. Construct a "super flat" by dividing the normalized Kurucz model into the normalized 
data (for both SP and S intervals). 
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6. Estimate the planet spectrum by subtracting the S interval spectrum from the SP 
interval spectrum. 

The four estimates of the exoplanet spectrum are then averaged to create the final spectrum. 
The errors are estimated by taking the average value of the differences in the estimate at 
each wavelength. 

To assess the magnitude of residual systematics, we made a comparison (see Fig. [TOl) 
between the spectra from each of the four independent AOR and nod combinations. In 
the central region of the IRS SLl instrument bandpass, the spectra are in relatively good 
agreement. However, this agreement becomes worse at either end of the instrument bandpass; 
this is especially true for wavelengths between 7.5 and 9 fim. The reason for this is that the 
assumption that G{X) is time invariant is only a first order approximation. Because we have 
normalized by the 12 fim fiux density, the effect of the small, uncorrected pointing errors 
within a nod is greatest at the band edges (see Fig. H]). To determine the best estimate of 
the differential spectrum, the four independent differential spectra are averaged together. 



3.4.2. Absolute Method: 

Here we describe how to apply the C correction to the source and calibrator data. 
Although we do not know a priori what the telescope pointing error is, we can determine the 
correct flux density for a given pointing error, y, using F = S/({y). From Eq. 4, we know 
that y = f{yo, niy, Ay, (py). Thus, our task is to identify the correct values for yQ, niy, Ay, (py. 
One way to do this is to require that the absolute spectrum be self-similar. We implemented 
this by constructing all unique combinations of the relation 

R(i,j)^F.(A.9)/F,(A.9,^giM (13) 

for the SP and S portions of each nod separately, where i and j are individual measure- 
ments in the time series. We then iteratively searched this space to determine the values of 
yo,my. Ay, (py (given in Tab. 2), which resulted in most closely approximating R(i,j) = 1. 
Fig. [7] shows the result of the application of the pointing offset correction, and the sec- 
ondary eclipse event is directly visible. Similarly, we applied the ( correction to t] Dor and 
determinded the pointing offset by requiring spectral self-similarity. As with the differential 
method, we evaluated the internal consistency of the pointing correction by comparing the 
spectra from both nods in both AORs. The agreement between the absolute spectra is ex- 
cellent (see Fig. [8]), and we now compare the differential and absolute spectra to assess the 
level of any residual systematics. 



-14- 



4. Discussion 

In this section, we compare tlie results of tlie differential and absolute methods we 
used to extract the planet spectrum. The assumption of wavelength-dependent stability 
used in the differential method is evaluated. We discuss our estimate of the eclipse depth 
and spectral features; we interpret these results in the context of recent models. We also 
discuss the significant d ifferences between our analysis methods and results and those of 



Richardson et al. I (120071 ). 



4.1. Comparison of differential and absolute methods 

A fundamental strength of our approach is the use of two semi-independent methods 
to demonstrate understanding and calibration of the dominant systematic errors. As Fig. 
[To] shows, the agreement between the differential and absolute planet spectrum estimation 
methods is excellent over most of the instrument passband. While agreeing within the 
errors, between 7.5 and 9 /im, the differential spectrum is systematically below the absolute 
spectrum. This is caused by small pointing errors occuring within a nod that are not removed 
by the differential method. Because the internal scatter of the absolute method is similar 
at all wavelengths, we consider the absolute spectra to be the best estimte of the planet 
spectrum. 

That the two spectral extraction procedures yield consistent results is encouraging and 
gives us a high degree of confidence that the calibration of systematic errors is successful 
within the error bars. Given that the differential method appears to make no specific cor- 
rection for pointing, one might wonder why the agreement with the absolute method is so 
good. The source of the agreement is that the normalization by the Kurucz model corrects 
for any chromatic error, pointing or otherwise, so long as the chromatic error changes in 
time are relatively small. Thus, normalization by the Kurucz model corrects the chromatic 
error produced by the largest pointing errors, which are static and occur during the initial 
peak-up and during the nod. Because the periodic pointing errors are relatively small, the 
change in the measured flux during a nod is, to first order, wavelength independent, and 
thus the spectral flat field is a good approximation for the flux correction due to the initial 
pointing error. Thus, the agreement between the differential and absolute spectral estima- 
tion methods supports the original assumption that the G{X) term is relatively (but not 
completely) constant during a nod. The increased size of the error bars in the differential 
method results from the periodic component of the pointing errors. 
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4.2. Eclipse Depth 



We have determined an average eclipse depth for the data by normalizing the absolute 
(pointing error corrected) flux density time series at each wavelength by the median values 
of the time series, F{X)/ < F{X) >. This is then averaged over wavelength to develop 
a broad-band light curve. The result of this can be seen in Fig. O the broad-band light 
curve clearly shows the eclipse and the transitions between ingress and egress. We can 
derive four independent estimates (one for each nod) of the broad-band eclipse depth, and 
these are consistent within the errors. After averaging the individual estimates, we flnd 
the average eclipse depth between 7.6 and 14 yum to be 0.00315±0. 000315. This minor 
restriction in wavelength was im plemented to exclude the channels with lower SNR. When 
compared to theoretical models (jSurrows et al. 2006 ). the measured eclipse depth suggests 
that substantial heat redistribution from the dayside to the nightside is occurring. This 
evidence of heat redistributi o n is s imila r to the interpre t ation given to observations of HD 



189733b bv lGrillmair et al. I (120071 ) and iKnutson et al. I fl2007h . 



4.3. Planet Spectrum 

We have determined the planet spectrum (see Fig. [TTj) and find it to range from about 
600 /iJy at 7.5 /im to about 200 /xJy at 15.2 /x. The SNR in the spectrum ranges from ~ 10 
at short wavelengths to ~ 2 at the longest wavelengths. To our knowledge, this is the first 
determination of the absolute spectrum of exoplanet emission. Our results for the spectral 
shape agree well with previous wo rk (see Fig. albeit wit h improved SNR; specifically, we 



confirm the marginal detection by iRichardson et al. I (120071 ) of a narrow feature near 7.7 fim. 



For most wavelengths, the planet spectrum is characterized by approximately featureless 
emission. However, between 7.5 and 8.5 /im, there is evidence for one broad (previously 
unreported) and one narrow (previously reported) spectral feature. 

Both the absolute spectrum and the contrast spectrum show evidence of a possible ~ 
0.5 /im-wide feature centered around 8.1 fim, with a significance of about 4 a. This broad 
feature represents a fiux deficit from the local trend and could be due to absorption. At 
the full spectral resolution, there is also a suggestion of a narrow feature around 7.7 /xm. 
This narrow feature candidate could be either in absorption (a deficit relative to the local 
trend in the 7.57 and 7.63 /im channels) or in emission (an excess relative to the local trend 
in the 7.69 and 7.75 /im channels). The shape of the broad feature causes us to favor the 
hypothesis of a narrow absorption feature in the 7.57 and 7.63 fim channels. However, 
movement of any one of these four spectral points (7.57, 7.63, 7.69, and 7.75 /im) by ~1.5 cr 
towards the local trend would convert this candidate feature into an outlier consistent with 
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a normal measurement distribution. The narrow feature candidate is sufficiently marginal 
that additional observations are required to con&m or rule out a spectral feature at this 
wavelength. 

The indication that the spectrum of HD 209458b contains one broad and one narrow 
feature between 7.5 and 8.5 /im is supported by the iRichardson et al. I (120071 ) measured 
spectrum. Indeed, the striking qualitative agreement (one broad and one narrow feature) 
between previous work and our results for the spectral mo dulation between 7.5 and 8 .5 iim is 
a strong indication that this modulation is real. Although IRichardson et al. I (120071 ) did not 
discuss the broad feature, it is present in their spectrum and we conffim their measurement. 
While we cannot totally exclude the possibility of some residual instrument systematic, it 
is highly significant that the shape of this spectral modulation is consistent using three 
independent methods conducted by two independent groups. Because of the repeatability of 
the result and the maturity of the exoplanet spectrum determination, the spectral modulation 
between 7.5 and 8.5 fim is likely real and may serve as a useful constraint on models for 
emission from HD 209458b. 



In the interpretation of the previous results for these data, IRichardson et al. I (120071 ) 



reported the detection of a broad emission feature centered at 9.65 fim, identified as a 
silicate feature, and a narrow emission feature centered at 7.78 fim. We find no evidence 
to support the identification of a 9.65 fim feature in our spectrum. Additional averaging 
and scrollin g median filtering does no t reveal any candidate feature with the characteristics 
claimed by 



Richardson et al. 



(120071 ). It is possible that the narrow feature identified by 
Richardson et al. I (l2007l ) corresponds to the 7.67 and 7.75 /xm channels in our analysis. If 
this is the case, the differenc e in wavelength is possibl y due to the non-standard wavelength 
calibration method used by IRichardson et al. I (120071 ) (see discussion below). However, we 
stress that the candidate absorption feature at 7.57 /xm is at least as likely as an emission 
feature at 7.69 fim 



4.4. Differences with Previous Work 



There are several significant dif ferences in our data calibr ation method and results when 
compared to the approach used by IRichardson et al. I (120071 ). Our approach explicitly cor- 
rects for the telescope pointing error and source leakage into the background; both of these 
effects are chromatic errors capable of introducting systematic errors in a spectrum. We also 
use two methods, one differential and one absolute, to extract the spectrum of the exoplanet, 
and we demonstrate good agreement between the methods. Unlike the previous work, we 
are able to explicitly measure the secondary eclipse depth from the IRS data. The improved 
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SNR and lower internal scatter in our spectrum allows a clear identification of the spectral 
modulation between 7.5 and 8.5 fim, and rules out the possibility of significant silicate emis- 
sion at 9.65 fim. Below, we explain so me of the important deta ils in the differences between 



our methods and results and those of iRichardson et al. I (120071 ) 



absolute spectrum (result): We have determined the spectrum of HD 209458b in 
Jy. To our knowledge, this is the first absolute determination of an exoplanet emission 
spectrum. 

eclipse depth (result): We explicitly determine the broad-band eclipse depth from 
the IRS data at high SNR (lOcr). This determines the eclipse depth in the IRS SLl 
instrument passband and avoids the uncertainty associated with incomplete matching 
of the IRS wavelength coverage to the 8 /im IRAC channel. 

spectral features (resu lt): We find no evidence for the silicate feature identified by 



Richardson et al. I (120071 ). There is the possibility of a narrow candidate feature at 
~7.7 /xm, but at the 1.5-cr level it is consistent with noise. In addition, the position 
of the 7.64 and 7.70 fim spectral points relative to the neighbors make this candidate 
feature as likely to be an absorption feature as an emission feature. 

• wavelength calibration (method): As part of the spectral extraction process, using 
SMART, we include the wavesamp.tbl table calibration file provided by the Spitzer 
Science Center. This approach implements an interpolation method to determine how 
fractions of a pixel contribute to a given wavelength. This approach accounts for the 
spectra tilt and curvature and provides Nyquist sampling of the spectra in the disper- 



sion direction. In contrast, the wavelength definition used by IRichardson et al. I (120071 ) 
is based on the bO_wavesamp_wave.fits file which, according to the IRS handbook, is 
for notional purposes only and should not be used for a scientific analysis. It is likely 
that relying on the bO_wavesamp_wave.fits file for the wav elength definition i s why 



the wavelength scales for the two AORs are different in the IRichardson et al. I (120071 ) 
analysis. 

background correction (method): Our background subtraction approach includes a 
correction for contamination from the source. This is a wavelength-dependent effect, 
which is of the order of the secondary eclipse depth. Failure to correct for source leakage 
in a normal background subtraction approach causes a wavelength -dependent error if 



the da ta in a nod are simply adjusted (the "multiplicative factor" for IRichardson et al. 



(120071 )) to make the time series continuous. 



- 18 - 



• pointing correction (method): Our method includes a specific correction for the 
pointing error, which corrects the static offset, periodic changes, and linear drift error 
terms in the telescope pointing. Uncorrected pointing errors that change with time 
introduce spectral errors. 

• spectral response function (method): Our determination of the spectral response 
function includes a correction for both the pointing error and the source contamination 
of the background. The spectral response function derivation is required for an absolute 
exoplanet spectrum. 

• error estimate (method): Our error bars are determined by the standard deviation 
in the mean of multiply determined quantities (e.g. the background corrected and 
pointing corrected ti me series) and by the roo t sum of squares for combined quantities. 



The error bars in the lRichardson et al. I (120071 ) analysis are determined by offsetting the 



time series by one time step, subtracting the original time series, and then determining 
the standard deviation in the mean of the resulting time series (in every spectral 
channel). This approach removes the effect of all systematic error with timescales 
longer than ~ 2 minutes and thus has the potential to underestimate the measurement 
uncertainty. 



Conclusions 



Our results for the spectrum of HD 209458b are consistent with a smooth, largely 
featureless spectrum ranging from about 600 /iJy at 7.5 /im to about 200 /xJy at 14 /xm. 
However, there is evidence of a spectral feature between 7.5 and 8.5 fim. We find evidence 
for a broad ~ 0.5 fim wide feature, centered at approximately 8.1 /im, that is possibly due to 
absorption. Near 7.7 fim we find a narrow feature candidate that could be either absorption 
or emission, depending on wavelength and local baseline trend assumptions; this candidate 
feature is only ~ 1.5 a from being cons i stent with noise. We find no evidence for the silicate 
feature reported in [Richardson et al. I (120071 ). The relatively smooth character of the HD 
209458b spectrum suggests the planet emission is dominated by purely thermal emission 
over most of the IRS SLl passband. However, the spectral modulation between 7.4 and 8.4 
fj,m is significant and suggests that t he dayside vertical te mperature profile of the planet 
atmosphere is not entirely isothermal (IFortney et al. II2006I ). 



We are able to make a direct measurement of the eclipse depth. Between 7.6 and 14.2 
/im we find an average ec l ipse d epth of 0.00315±0. 000315; when compared to planet emission 
models such as iBurrows I (120071 ). the measured eclipse depth is suggestive of substantial heat 
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redistribution between the ni ghtside and dayside. Similar conclusions ha ve been drawn for 



observations of HD 189733b flGrillmair et al. 1120071 : iKnutson et al. 1120071 ) 



The methods we have developed for calibration of the background and pointing errors 
represent a significant improvement in the state of the art for IRS calibrations on bright 
objects. Using a simple pointing model and requiring self-consistency of the spectrum for 
the "star+planet" and "star" portions of the time series, we are able to optimally recover 
the spectrum of HD 209458b. By applying our calibration of (i) source contribution to the 
background and {ii) pointing errors to the definition of the spectral response function, we 
have achieved an absolute flux density calibration approaching 0.1 %. This implies that our 
calibration method is suitable for spectroscopy of emission from the nightside of exoplanets 
and would significantly increase the SNR for IRS spectra of relatively bright point sources. 



We thank the original PI team for the proposal and preparation of the AORs required 
to obtain these data. We appreciate the comments of the anonymous referee who encouraged 
us to fully describe our calibration process and extend our calibration method to the entire 
IRS SLl passband; these suggestions led to the detection of spectral modulation that was 
outside our original spectral passband. We also thank Drake Deming for several helpful 
conversations regarding the reduction of secondary eclipse data. We thank John Bayard for 
several helpful discussions concerning Spitzer pointing errors and Sara Seager for discussion 
regarding the possible role of clouds in exoplanet atmospheres. 
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Fig. 1.— The 9 /xm flux time series for AOR 14817792 and AOR 14818048. In botli AORs, 
a nod was executed near tlie center of tlie secondary eclipse period. The measured flux is 
affected by systematic errors caused by (i) pointing, (ii) background subtraction, and (iii) 
latent charge accumulation; these effects combine to cause changes in the measured flux that 
are ~10x larger than the secondary eclipse signature. The pointing error causes the offset 
between the flux density in the nod 1 data in both AORs, and it causes most of the offset 
in flux density between the nod 1 and nod 2 data for a given AOR; the pointing error also 
causes the linear and periodic changes in the measured flux. The background subtraction 
error causes a portion of the offset between the nod 1 and nod 2 data for a given AOR. The 
latent charge accumulation error causes the initial rapid increase in measured flux; this 
effect is larger in AOR 14817792 due to the array history. 
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Fig. 2. — A schematic representation of the steps in our data reduction procedure; a detailed 
discussion of each stage of this process is contained in the text. The strength of this approach 
is that it provides two semi-independent methods of estimating the planet spectrum. Because 
the two methods interact with systematic errors differently, the comparison of the results 
provides a robust method for demonstrating control of these errors. In the final step of the 
data reduction, the estimates of the planet spectra from the two methods are compared. 
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Fig. 3. — Top: Source position data in the slit spatial axis, together with an elliptical 
pointing model for AOR 14818048. Both periodic motion and a linear drift are clearly 
identifiable in the source position data. The data are fit to determine the amplitude, period, 
eccentricity, and linear drift terms. After determining the best fit pointing model using the 
data in the spatial axis, we investigated the effect of the cross-slit motion predicted by this 
model. Bottom: The normalized cross-slit projected motion from the elliptical pointing 
model as a function of time; with the exception of the initial few points, the excellent 
agreement between the model and the data is evident. This strongly suggests that time- 
variable pointing errors (periodic changes and linear drifts) are primarily responsible for 
changing the measured flux. 
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Fig. 4. — Top: The change in measured flux as a function of source position in the cross-sht 
direction of the IRS entrance sht. This cahbration was developed from sht scans of IRS 
calibrators (see text). The measured flux does not peak at zero offset because of an asym- 
metry in the wings of the PSF. Given the 0.4 arcsecond radius of the IRS "high precision" 
peakup mode error circle, the uncertainty associated with the standard IRS high-precision 
peakup results in a ~ 7 % uncertainty in the measured flux. The measured flux error due to 
the cross-slit pointing error is, to flrst order, common mode in wavelength. Bottom: The 
pointing offset correction as a function of wavelength normalized to 12 //m. This plot shows 
that there is a spectral error of ~ 1.5 % across the instrument passband due to the pointing 
offset effect for the standard high-precision peakup. The horizontal scale is set to include 
the range of offsets present in the observations of HD 209458. 
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Fig. 5. — Top: The background source contribution function (BSCF) is different for the 
two nods because of an asymmetry in the telescope PSF (see text for description of the 
derivation). Multiple stars and epochs were used to determine both the correction and the 
uncertainty. Bottom: The open circles (nodi background) and triangles (nod2 background) 
are the estimated background prior to correction for the source contribution at 12 /im in 
AOR 14818048. The filled circles and triangles show the result of applying the correction 
for source leakage. 



- 26 - 



I 




Wavelength (microns) 

Fig. 6. — To determine the level of internal consistency for the differential method, we have 
plotted the four independent estimates of the planet's spectrum. The internal agreement is 
good over most of the instrument bandpass. However, this internal agreement is substantially 
worse between 7.5 and 9 /xm; this is due to small, uncorrected pointing effects that are 
greatest at the edges of the passband because we have normalized by the 12 /xm flux density. 
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Fig. 7. — The 12 fim flux time series for AOR 14817792 and 14818048 as a function of orbital 
phase. The background-subtracted data are shown in open circles and triangles. The large 
discontinuity halfway through each data set corresponds to a nod; periodic changes and drifts 
also affect the time series measurements. The changes in measured flux are primarily due to 
telescope pointing errors, which cause variable vignetting by the spectrometer entrance slit of 
the wings of the source point spread function. The filled points show the pointing-corrected 
flux. The signature of the secondary eclipse is visible as an interval with systematically lower 
flux density. The vertical dashed lines indicate the region of ingress/egress where the planet 
is partially obscured by the star. A horizontal line has been added to aid in identifying 
the flux density decrease during secondary eclipse. Note that the effect of latent charge 
accumulation can be seen in the apparently increasing flux in the initial ~ 20 points for each 
AOR (near orbital phase 0.464). 
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Fig. 8. — Four independent estimates (nodi and nod2 for both AOR 14817792 and AOR 
14818048) of the planet spectrum, determined using the absolute method, are plotted to- 
gether as an internal consistency check. The internal agreement of the absolute method is 
excellent over the entire instrument passband. 
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Fig. 9. — The broad-band normalized light curves average between 7.6 and 14.2 /xm. The 
filled points are AOR 14817792, and the open points are AOR 14818048. The eclipse is 
clearly evident, as are the ingress/egress transitions. The rapid increase in the points at 
the beginning of each AOR shows the accumulation (and subsequent stabilization) of latent 
charge. We have omitted the data affected by latent charge (indicated by the shaded region) 
from our analysis, and we find the average echpse depth is 0.00315 ±0.000315. 
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Fig. 10. — The mid- infrared spectrum of HD 209458b derived using the absolute method 
(circles) and the wavelength differential method (triangles). For both cases, the spectrum 
has been normalized to the 12 /xm flux and scaled so that the average flux for both plots is 
identical. While the agreement is excellent over most of the passband, the effect of small, 
uncorrected pointing errors causes the differential spectrum to lie slightly below the absolute 
spectrum at some wavelenghts between 7.5 and 9 /xm. 



-31 - 




8 10 12 



Wavelength (microns) 

Fig. 11. — The spectrum of HD 209458b is mostly smooth but contains some modulation 
between 7.5 and 8.5 fim. There is the suggestion of a spectral feature at about 14.5 fim; 
however, the SNR in this region of the spectrum is low, and the location of the spectral 
baseline is uncertain. 
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Fig. 12. — The contrast spectrum for H P 209458b as deter r ained by our calibration method 
(black circles) and the method used by iRichardson et al. I (120071 ) (red triangles). Over the 
range of wavelengths where the spectra can be directly compared, the internal scatter of 
our spectrum is 2.3 times smaller than the previous result. The methods are in qualitative 
agreement, and both show evidence for modulation of the spectrum between 7.5 and 8.5 /xm. 
This modulation includes a broad feature seen as a local minima centered at approximately 
8.1 fim, and a narrow feature that could be either a local minima at approximately 7.57 fim 
or a local maxima centered at approximately 7.63 fim. 
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Fig. 13. — The contrast spectrum for HD 209458b plotted together with a model (Burrows 
et al. 2006) for the exoplanet emission. The planet emission is most consistent with models 
containing significant heat redistribution. The departure of the measurement from the model 
in the 8 /im region of the spectrum is significant and may place constraints on models of the 
planet emission. We have omitted wavelengths beyond 13 fim from the contrast spectrum 
because these data have lower SNR. 
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Table 1: The parameters determined by fitting an elliptical pointing error to the source 
position along the spatial axis of the spectrometer slit. Aq is the amplitude of the angular 
acceleration, uj is the angular frequency, 0e is the angular phase, Xq is the constant offset in 
the X position, rux is the linear change of position with time in the spatial axis, is the 
amplitude of the x periodic function, and 0^:, is the phase in the spatial axis. 



Spectral Axis Fit Parameters 



AOR/nod 


yo 


my 


Ay 


4>y 




[pixels] 


[pixels/nod] 


[pixels] 


[phase] 


AOR 14817792/nodl 


0.2149 


-0.0242 


0.0145 


0.8474 


AOR 14817792/nod2 


0.3003 


-0.0276 


u 


a 


AOR 14818048/nodl 


0.1175 


-0.0202 


u 


a 


AOR 14818048/nod2 


0.2023 


-0.0403 


u 
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Table 2: The parameters for the telescope pointing model in the cross-slit (dispersion) axis: 
initial offset, t/o! linear drift, my\ periodic amplitude, Ay\ and phase, (py. These values were 
determined by requiring self-similarity of the spectrum during the "star-|-planet" and "star 
only" intervals of the spectral time series. 
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Planet/Star Contrast Spectrum 



wavelength (/i) 


contrast 


error 


7.67 


0.0031 


0.00025 


8.15 


0.0027 


0.00019 


8.63 


0.0032 


0.00020 


9.12 


0.0033 


0.00022 


9.60 


0.0032 


0.00022 


10.09 


0.0031 


0.00021 


10.57 


0.0033 


0.00023 


11.05 


0.0034 


0.00026 


11.54 


0.0030 


0.00028 


12.02 


0.0033 


0.00029 


12.51 


0.0032 


0.00029 


12.99 


0.0033 


0.00033 


13.47 


0.0025 


0.00040 


13.96 


0.0029 


0.00046 


14.44 


0.0040 


0.00049 



Table 3: Data and l-cx error for the contrast spectrum of HD 209458b shown in the lower 
panel of Fig. 13. 



